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Abstract 

We discuss short-baseline and very-short-baseline v e disappearance at a neutrino factory. 
We take into account geometric effects, such as from averaging over the decay straights, 
and the uncertainties of the cross sections. We follow an approach similar to reactor exper- 
iments with two detectors: we use two sets of near detectors at different distances to cancel 
systematics. We demonstrate that such a setup is very robust with respect to systematics, 
and can have excellent sensitivities to the effective mixing angle and squared-mass splitting. 
In addition, we allow for CPT invariance violation, which can be tested (depending on the 
parameters) up to a 0.1% level. 
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1 Introduction 



Neutrino oscillation experiments have shown that neutrinos are massive particles with at 
least two squared-mass differences: Amg OL ~ 8 x 10~ 5 eV 2 , measured in solar and very-long- 
baseline reactor neutrino experiments, and Am\ TM ~ 2 x 10~ 3 eV 2 , measured in atmospheric 
and long-baseline neutrino experiments (see Refs. [1-10]). These two Am 2 's are perfectly 
accommodated in the framework of three-neutrino mixing, where there are two independent 
squared-mass differences. However, there are experimental anomalies which may indicate 
the existence of Short-BaseLine (SBL) or Very-Short-BaseLine (VSBL) oscillations gener- 
ated by a third Am 2 which is much larger than the other two: Amg BL > 10 _1 eV 2 or 
Amy SBL > 10 eV 2 . Among these anomalies, the most well-known is the LSND signal in 
favor of SBL — > v e oscillations [11], which has not been confirmed by other experiments 
and is currently disfavored by the negative results of KARMEN [12] and MiniBooNE [13]. 
Less well-known are the Gallium radioactive source experiments anomaly [14] and the Mini- 
BooNE low-energy anomaly [13], which could be explained by SBL [15,16] or VSBL [17,18] 
v e disappearance. 

The existence of a third Am 2 requires the existence of at least a fourth massive neutrino 
which corresponds, in the flavor basis, to the existence of a sterile neutrino fermion 
which is a singlet under the Standard Model symmetries. Hence it is electrically neutral and 
does not take part in weak interactions. If the three active neutrinos z/ e , u^, and v T are mixed 
with the sterile neutrino, neutrino oscillation experiments can observe the disappearance of 
active neutrinos into v s . 

In light of the above-mentioned anomalies, it is interesting to investigate the possibility of 
(V)SBL v e disappearance with future high-precision experiments. In general, it is important 
to investigate the possibility of v e disappearance generated by a Am 2 different from Am| OL 
and Am 2 ATM in order to constrain schemes with mixing of four (see Refs. [1,4,6,8]) or 
more [19,20] massive neutrinos. These schemes have been studied mostly in connection with 
the LSND anomaly, but the latest global fits of the experimental data, including the LSND 
signal, are not good [4,8]. However, the schemes with mixing of more than three neutrinos 
may be realized in nature independently of the LSND signal. Hence, it is important to 
investigate the phenomenology of sterile neutrinos with an open mind, not only through 
neutrino oscillations [21-29], but also by studying their effects in astrophysics [30-35] and 
cosmology [36-38]. 

If there is (V)SBL electron neutrino disappearance, it must be mainly into sterile neutrinos, 
because the mixing of the three active neutrinos with the fourth massive neutrino must 
be small in order to fit the data on v e — > v^^ oscillations generated by Am| OL and the 
data on — > v T oscillations generated by Am 2 ATM . In the 3+1 four- neutrino schemes (see 
Refs. [1,4,6,8]) with Am^ SBL = lAm 2 ^ > Am 2 ATM = lAmf^ » Amg OL = (Amiy, 
where Am 2 k - = m 2 — m 2 , the mixing matrix U must be such that \U e ^\, {U^l, \U T ±\ <C 1 and 
\U S i\ ~ 1. Therefore, the amplitudes of the (V)SBL oscillation channels, A a p = 4\U a ^\ 2 \U/34\ 2 
for a ^ /3, are such that A ab A as for a,b = e, /i, r. 

In this paper we study the sensitivity of neutrino factory experiments to (V)SBL v e and v e 
disappearance, which in practice has been investigated so far mainly through SBL reactor 
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neutrino experiments [y e disappearance). 

We will first study, in Section HI (V)SBL v e and v e disappearance at a neutrino factory- 
assuming exact CPT symmetry, which implies P ee = P g g (see Ref. [7]), considering the 
simplest case of effective two-neutrino mixing with 

P ee = P m = 1 - sin 2 (2£) sin 2 , (1) 

where, from now on, Am 2 = Am 2 v ^ SBL . This is the case of four-neutrino mixing schemes 
with Am 2 = lAmlJ ^> Am^ TM = [Am 2 ]] 3> Amg OL = jAm^]. In the 3+1 schemes, 
the amplitude of the oscillations is related to the U e 4 element of the mixing matrix by 
sin 2 (2fl) = A\U e4 \ 2 (l - \U e4 \ 2 ) (see Refs. [1,4,6,8]). 

The CPT symmetry is widely believed to be exact, because it is a fundamental symmetry 
of local relativistic Quantum Field Theory (see Ref. [39]). However, in recent years studies 
of extensions of the Standard Model have shown that it is possible to have violations of 
the Lorentz and CPT symmetries (see Refs. [40-42]) and several phenomenological studies 
of neutrino oscillations with different masses and mixing for neutrinos and antineutrinos 
appeared in the literature [43-54]. We will consider this scenario in the simplest case of 
effective two-neutrino mixing with 

P ee = 1 - sin 2 (2^) sin 2 (^f^) , (2) 

P ss = 1- S in 2 (2^)sin 2 (^^) . (3) 

This kind of CPT violation in a four-neutrino mixing scheme could reconcile the LSND signal 
with the other neutrino oscillation data [50] and/or could explain the Gallium radioactive 
source experiments anomaly and the MiniBooNE low-energy anomaly together with the 
absence of v e disappearance in reactor neutrino experiments [18]. Let us emphasize that the 
reconciliation of the LSND anomaly with the results of other neutrino oscillation experiments 
is not possible in three-neutrino mixing schemes even if CPT violation is allowed [8,52]. 

Another hint in favor of a possible CPT violation comes from the recent measurement of 
and disappearance in the MINOS experiment [55], which indicate different best-fit values 
of the oscillation parameters of and v^. Am^ INOS ~ 2 x 10~ 2 eV 2 and sin 2 #minos — 0.6 
for v^s, whereas A?71m INOS ~ 2.4 x 10~ 3 eV 2 and sin 2 Aminos — 1 f° r *V S - The best-fit val- 
ues and allowed region of the oscillation parameters are in agreement with atmospheric 
v ti v T oscillations. Since the 90% C.L. allowed region of the oscillation parameters 
has a marginal overlap with the much smaller 90% C.L. of the oscillation parameters 
(see the figure in page 11 of Ref. [55]), the MINOS hint in favor of CPT violation is rather 
speculative. Nevertheless, it is interesting to notice that a global separate analysis of neu- 
trino and antineutrino data in the framework of three-neutrino mixing with CPT violation 
leads to different best-fit values of the oscillation parameters of neutrinos and antineutrinos 
with Am 2 ATM ~ AmM INOS and sin 2 $atm - sin 2 Aminos, whereas Am 2 ATM ~ AmJj INOS and 
sin 2 #atm — sin 2 Aminos [56]. However, in this paper we do not consider the MINOS hint in 
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favor of CPT violation. We concentrate our study on possible CPT violations in (V)SBL 
v e and v e disappearance due to squared-mass differences larger than about 0.1 eV 2 . 

Besides those in Eqs. (j2J) and ([3]), it is possible to consider other, more complicated, ex- 
pressions for P ee and P ss , with additional energy-dependent terms in the oscillation phases 
which could be generated by modified dispersion relations that are different for neutrinos 
and antineutrinos (see, for example, Refs. [57-62]). However, the introduction of more 
unknown parameters would make the analysis too cumbersome, without much additional 
information on the potentiality of a neutrino factory experiment to test CPT invariance. 
In fact, it is plausible that the additional energy-dependent terms in the oscillation phases 
generate spectral distortions which would make the identification of new physics even easier 
than in the simplest case that we consider. 

In order to test CPT invariance (or small deviations from it) explicitly, it is convenient to 
define the averaged neutrino oscillation parameters 

9 = i {9 V + 9 9 ) , Am 2 ^ 1 - (Ami + Am 2 ) , (4) 
together with the CPT asymmetries 

_9 V -9 V _ Ami ~ / K x 

° CPT = ^T^' mCPT = Ami + Am? ' (5) 

which are constrained in the range between —1 and 1. Then we have 

9 V = (l + a CPT )9, (6) 

9 D = (l-ac PT )0, (7) 

Ami = (l+m CPT )Am 2 , (8) 

Ami = (! - ™cpt) Am 2 . (9) 

The limit of CPT invariance (Eq. ([!])) corresponds to ocpt = wcpt = 0. In Section [5] 
we discuss the potentiality of neutrino factory experiments to discover acpT ^ and/or 
m C pT 0. 

The plan of the paper is: in Section [2] we define an "ideal detector" for the measurement 
of (V)SBL v e and v e disappearance at a neutrino factory, and we describe our treatment of 
geometric effects; in Section [3] we discuss the requirements for systematics; in Section 0] we 
discuss the sensitivity to (V)SBL v e and v e disappearance assuming CPT invariance, with 
the survival probability in Eq. (JT]); in Section [5] we discuss the sensitivity to CPT violation 
considering the survival probabilities in Eqs. ([2]) and ([3]); conclusions are presented in the 
final Section [6j 

2 Ideal detector and geometric effects 

Our neutrino factory geometry is based on the International Design Study for the Neutrino 
Factory (IDS-NF) baseline setup [63], with the geometry illustrated in Fig. HJ We consider 
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Figure 1: Geometry of the decay ring (not to scale). Two possible detector locations are shown at 
d = 50 m and d = 2 000 m, where d is the distance to the end of the decay straight. The baseline L is the 
distance between production point and detector. 

2.5 x 10 20 useful muon decays per polarity and year, with muon energy = 25 GeV. For 
the total running time, we consider ten years. 

In order to test SBL v e disappearance, we add detectors in front of the decay straights as 
illustrated in Fig. HJ Here "near" and "far detectors" refer to SBL v e disappearance only, 
whereas the detectors for standard oscillations are much farther away and not relevant for 
our problem. The straight sections are anticipated to be about s = 600 m long. The distance 
d is the distance between the end of the decay straight and the near detector. The baseline 
L is the distance between production point and near detector, i.e., d < L < d + s. Since 
the n + and \x~ are assumed to circulate in different directions in the ring, we need pairs of 
detectors in front of the straights because we want to test CPT invarianceQ 

Since there are no specifications for near detectors at a neutrino factory yet (see Ref. [64] 
for a generic discussion), we turn the argument around and formulate the requirements 
for the detectors for this measurement. Our detectors are assumed to measure the total 
charged current rates with a 100% detection efficiency; a lower efficiency will simply lead 
to a re-scaling of statistics and can be easily compensated by a larger detector mass. The 
energy threshold is chosen to be 500 MeV, similar to a Totally Active Scintillator Detector 
or an iron calorimeter, and the energy resolution is taken as 



with e = 0.55 GeV and Eq = 1 GeV, which is a conservative estimate for a magnetized iron 
calorimeter [63] . Similarly, we assume that the neutral current level can be controlled at the 
level of 10 -3 from all neutrinos in the beam (see, e.g., Refs. [65,66] in the context of a low 
energy neutrino factory). However, we have tested that the results do not strongly depend 
on these three quantities. We require an excellent flavor identification (at the level of 10 -3 
for the misidentification, as we will see later). Charge identification is also desirable in 

1 Without CPT invariance test, detectors in front of one straight are sufficient. The detectors in front of 
the other straight only increase statistics then. 




(10) 
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order to reduce the contamination of the v e (or v e ) signal by v e (or u e ) generated by possible 
(V)SBL — > v e (or — > u e ) oscillations. However, we do not consider the backgrounds 
from charge misidentification explicitly!! For the binning, we use 17 bins between 0.5 and 
25 GeV with a bin size of 0.5 GeV (1 bin) - 1 GeV (9 bins) - 2 GeV (5 bins) - 2.5 GeV (2 
bins). As the main obstacles for the physics potential, we have identified the extension of 
the decay straights and the impact of systematics. We discuss the first issue below, and the 
second issue in the next section. Thereby, we define our "ideal detectors" as detectors with 
the above properties, but no backgrounds and systematics. 

Our geometric treatment of the near detectors is based on Ref. [67], which discusses the 
flux at near detectors in detail. Here we start from the differential event rate from a point 
source dNps/dE without oscillations. Taking into account the extension of the straight and 
the geometry of the detector, the averaged differential event rate is given b}{§ 

rfiVavg 1 d fdN 1 TdNpsi^E) 

-dL = - I — e(L, E) P ee (L, E)dL . (11) 



dE s J dE s J dE 

d d 

Here e(L,E) = A e g/Ai) et parameterizes the integration over the detector geometry for a 
fixed baseline L and given energy E (v4rj e t is the surface area of the detector and A c q is 
the effective surface area which takes into account the angular dependence of the neutrino 
flux). Since dNpg/dE oc 1/L 2 , we can re- write this as 

rfiV avg _ diVps(L e ff, E) Lis d Te(L,E) , r dN PS (L c{i ,E) 



dE dE s 

d 



l^ Pte{L , E)dL = ^^lp (E) , (12) 



with the average efficiency ratio times probability!^ 



d+s 

L 2 cS fe(L,E) 



P (E) = -2* I P ee (L, E)dL , (13) 



and the effective baseline 



^/d(d+s), (14) 

such that P(E) = 1 for e(L,E) = P ee (L,E) = 1. We assume e(L,E) = 1 (far distance 
approximation), which, to a good approximation, is satisfied for ND4 of Ref. [67] (see Fig. 4 
therein) for d > 50m. This detector is very small (200 kg), however, with a sufficient event 
rate. At a neutrino factory, the active volume of near detectors are probably going to be 
rather small, because high granularity and good track reconstruction will be more important 



2 The level of contamination depends on the oscillation model. Even for large mixing angles driving these 
oscillations of the potential background, a charge misidentification level of about 10~ 3 would be sufficient. 

3 Notc that as a peculiarity compared to Ref. [67], dNps/dE uses the unoscillated event rate, because 
the oscillation probability has to be integrated over. 

4 Note that Eq. t|12J) implies that in GLoBES a point source spectrum at the effective baseline L c s can 
be used, which has to be corrected by Eq. (Tl3|) . We perform Eq. t|13|) directly in the probability engine. 



5 




sin 26 



Figure 2: Exclusion limit for several near detector distances d and our ideal near detectors (CPT 
invariance assumed; 90% CL, 2 d.o.f.; two near detectors in front of straights). The dashed curves illustrate 
the effect of including the averaging over the decay straight, whereas the solid curves are without this 
averaging. The fiducial detector masses are fixed to 200 kg. Note that there is no systematics included in 
this figure. 



than the active volume size [64]. Our "ideal" test detectors therefore have 200 kg fiducial 
volume at very short distances. One can, for longer baselines, up-scale the detector mass as 



d x (d + 600 m) n n 

m Det ~ ^O^t 

50 m x 650 m 



(15) 



without strong geometric effects from the effective area of the detector (i.e., one still operates 
in the far distance limit). However, one may choose a different technology for these larger 
detectors. 

For our simulation, we use the GLoBES software [68,69]. We define the exclusion limit as 
a function of sin 2 29 and Am 2 as the excluded region obtained in a \ 2 analysis assuming 
a vanishing true value of 9 (i.e. no oscillations). In Fig. [21 we show this exclusion limit 
for several near detector distances including the effects of averaging over the decay straight 
(dashed curves) and without averaging (solid curves). This figure is based on our "ideal" 
detectors without taking into account systematics yet. Obviously, the optimal detector 
locations depend on the region of sensitivity of Am 2 which is of interest: the smaller Am 2 , 
the longer the baseline. For instance, for Am 2 ~ 1 eV 2 , best sensitivity is obtained for 
d ~ 20 km, whereas for Am 2 ~ 100 eV 2 , a distance of the order d = 100 m is optimal. For 
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short distances d up to a few hundred meters, there is clearly an effect of the averaging 
over the decay straight. However, note that because of the 1/L 2 weighting in Eq. (fl2|) . the 
effect becomes negligible for d > 1 km. Compared to a classical beam dump experiment, one 
cannot get arbitrarily close to the source without loosing information. In the next section, 
we will discuss the requirements for systematics. 

We have also tested a low energy neutrino factory for this measurement, with similar success. 
However, in the absence of official numbers for the storage ring geometry and systematics, 
we will not discuss it in greater detail. In addition, note that the absolute performance 
is not a priori better than for a higher energy neutrino factory. For instance, assume 
that the distance d is fixed for geometry reasons. Then the oscillation effect is, to a first 
approximation, proportional to 1/E 2 (with E the peak energy of the spectrum), but the 
statistics roughly increases as E 3 (E 2 from the beam collimation and E from the cross 
sections), which means that the net effect is proportional to E^. We observed this behavior 
in our simulation. 

3 Requirements for systematics 

As far as systematics is concerned, it is well known from reactor experiments, such as 
Double Chooz [70] and Daya Bay [71] , that electron neutrino disappearance is most affected 
by the signal normalization uncertainty (see, e.g., Refs. [72,73]). We expect the same for 
our measurement. However, compared to reactor experiments, our signal normalization 
error does not mainly come from the knowledge on the flux, which we may know to the 
level of 0.1% using various mean monitoring devices [64], but from the knowledge of the 
cross sections. Because our neutrino energies span the cross section regimes from quasi- 
elastic scattering, over resonant pion production, to deep inelastic scattering, it is not a 
priori simple to estimate the accuracy of the cross sections knowledge at the time of the 
measurement. For reactor experiments, on the other hand, the inverse beta decay cross 
sections are well known. Note that Ref. [29] also uses this well-understood detection reaction 
for a low-gamma beta beam, whereas we will use a completely orthogonal strategy. 

Let us first of all illustrate what the main requirements for systematics are. As indicated 
above, we have tested in Fig. [3] the impact of a signal normalization error and an additional 
tilt error (tilting the shape of the spectrum). Although the errors are assumed to be rather 
optimistic (2.5%), there is a significant impact on the sensitivities at all baselines, as we 
expected. Off the oscillation maxima, as visible in the right panel at large values of Am 2 
where P ee ~ 0.5 sin 2 29, the signal normalization error ONorm directly limits the sensitivity 
to sin 2 29 ~ 2 V2.3 o"Norm — 0.076 at la (2.3 is the A% 2 corresponding to la for two degrees 
of freedom). The tilt error tilts the spectrum linearly, and is a first order approximation 
for a shape error. It is especially important where the spectral information leads to a good 
sensitivity, in particular, for the shorter baselines (left panel). However, note that this (lin- 
ear) tilt error cannot fully take into account the uncertainties in the cross sections, because 
the actual deviation may be non-linear. We have also tested the impact of backgrounds, 
energy resolution and energy threshold. The most important of these three systematics is 
the background, where the sensitivity is basically limited by the product of background 
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Figure 3: The effect of a different (hypothetical) systematical errors: A signal normalization error of 
2.5% and an additional spectral tilt error of 2.5% has been applied to the exclusion limit for two different 
detector distances d (90% CL, 2 d.o.f.). The dashed curves refer to our ideal detectors, the solid curves 
include systematics. Here the fiducial mass is fixed to 200 kg, the effect of averaging over the decay straight 
is taken into account. Here CPT invariance is assumed. 

level and background uncertainty. Even for large uncertainties of the background, such as 
20%, this product limits the sensitivity to about 0.001 x 20% ~ 10~ 4 , which is beyond our 
expectations in the presence of a normalization uncertainty. 

In summary, the signal normalization and shape have to be either very well known, or very 
well measured. The first requirement means that one needs very refined theoretical models 
for the cross sections, the second possibility means that one needs to measure the cross 
sections very well. We follow the second approach by considering a setup with two sets of 
detectors (c/., Fig. [I]): 

1. Near detectors at d = 50 m with mDot = 200 kg. 

2. Far detectors at d = 2 000 m with m-Det = 32 1. 

The signal measured with the near detectors fixes the normalization and shape of the un- 
oscillated signal (for small enough Am 2 ). The far detectors are up-scaled versions of the 
near detectors following Eq. (]15[) . which means that geometric effects are almost negligible. 
The near detectors have optimal sensitivity at a few hundred eV 2 (VSBL), whereas the far 
detectors have optimal sensitivity at a few eV 2 (SBL). Note that longer baselines may be 
even better for the far detectors, but then the depth difference between storage ring and 
detectors may become unrealistically large. On the other hand, for distances much shorter 
than 2 km, one significantly looses sensitivity for small Am 2 . 
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For systematics, we adopt the most conservative point of view, i.e., we assume that we 
hardly know anything about the cross sections, neither the normalization nor the shape, 
but that the cross sections are fully correlated among all detectors measuring the cross 
sections. Such an error is often called "shape error" and is uncorrelated among the bins. In 
summary, we include the following systematical errors similar to the reactor experiments in 
Ref . [73] , and we have tested their impact (we have switched off systematical errors to test 
their impact): 

Shape errors uncorrelated among bins and v-v, but fully correlated among the detectors. 
These errors include cross section errors, scintillator or detector material properties, 
etc.. In addition, flux errors can be included here (the detectors only measure the 
product of flux and cross section for the disappearance channel). We estimate this 
error to be 10%. However, even a larger error does not matter if both near and far 
detectors are present, but only errors considerably smaller than 10~ 3 improve the 
result significantly (which is absolutely unrealistic for this type of systematics). 

Normalization errors uncorrelated between the near and far detectors. These relative 
normalization errors come from the knowledge on fiducial mass, detector normaliza- 
tion, and analysis cuts (uncorrelated between the detectors) . They are typically small 
if similar detectors are used. For reactor experiments (Double Chooz [70]), this error 
is about 0.6%, which we use as an estimate. We have tested that there is little depen- 
dence on this error unless it can be reduced to the level of 10~ 4 (then there is a small 
improvement), if the other systematics is present. 

Energy calibration errors uncorrelated between the near and far detectors of the order 
0.5% are used (similar to the reactor experiments). As we have tested, they are of 
secondary importance if all the other systematics is present. 

Backgrounds at the level of 10~ 3 from neutral current events etc. are assumed, known to 
the level of 20% (somewhat conservative estimate from a neutrino factory). If all the 
other errors are present, backgrounds hardly matter. 

The effect of electron neutrino disappearance on the event rates of the individual bins is 
illustrated in Fig. H] for the near (left) and far (right) detector for several values of Am 2 . 
For relatively small Am 2 ~ leV 2 (diamond curves), the near-far combination will perform 
similar to the reactor experiments with two detectors, where the near detector measures 
the shape and the far detector the oscillation effect. For Am 2 3> 1000 eV 2 (c/., triangle 
curves for comparison), the oscillations average out in both detectors, and sin 2 29 can only 
be constrained to the level of the shape errors (whereas Am 2 cannot be measured). For 
Am 2 ~ 100 eV 2 (box curves), the oscillation effect will mainly take place in the near detector, 
whereas the far detector measures the shape (after averaging). For Am 2 ~ 10 eV 2 (star 
curves), the situation is most complicated: there are oscillation effects in both detectors, 
which can lead to intricate parameter correlations. 
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Figure 4: Relative effect on the binned (neutrino) event rates for several values of Am 2 , and sin 2 29 = 0.1, 
in the near (left) and far (right) detectors. For each energy bin we plotted (R — Rq)/Rq where R and Rq 
are the expected rates with and without oscillations. 

4 Results for CPT invariance 

All results presented in this section are based on our two-baseline setup without refined 
systematics treatment, assuming CPT invariance, i.e., the equal electron neutrino and an- 
tineutrino survival probabilities in Eq. (TO . 

Fig. [5] shows the performance of our near- far model (thick curve), where the effect of using 
only one set of detectors (near or far) is shown separately as thin curves. If only one set of 
detectors is used, the result will be limited by the 10% shape errors, i.e., it depends on the 
assumptions used. However, if the two sets of detectors are used, the impact of systematics 
cancels and the result is very robust with respect to the assumptions. From the above 
discussion, it should be clear that the results in this case do not depend very much on the 
actual numbers for the systematical errors. Nevertheless there is a considerable deviation 
from the no-systematics case (dashed curve). The improvement towards this hypothetical 
sensitivity requires a very good understanding of the cross sections at the level of the sin 2 29 
sensitivity. We have also checked that the performance cannot even be significantly improved 
with considerably larger detectors, because of the systematics limitation (even without the 
geometric effect of the beam included). 

Figure [5] shows that the sensitivity of a neutrino factory experiment to (V)SBL v e disap- 
pearance represents a dramatic improvement with respect to the sensitivity of reactor ex- 
periments, which is at the level of sin 2 29 ~ 10" 1 at large values of Am 2 (c/., thin gray/cyan 
curve). Moreover, the neutrino factory measurement with the near-far detector setup dis- 
cussed in Section [3] is model-independent, whereas reactor measurements of P ss depend on 
the calculated flux of u e 's produced in a reactor. Reactor neutrino experiments cannot take 
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Figure 5: Exclusion limit in the sin 2 29- Am 2 plane for our default configuration including systematics 
(thick solid curve, 90% CL, 2 d.o.f.). The thick dashed curve refers to our ideal detectors (no systematics), 
with near (ND) and far (FD) detectors combined. The thin solid curves illustrate the results for the near 
(50 m) and far (2 km) detectors if operated separately, but with full systematics. The effects of averaging 
over the decay straights are taken into account. The thin dashed curve corresponds to the default beta 
beam setup from Ref. [29] for comparison. The thin gray/cyan curve is the current limit from Bugey [74] 
+ Chooz [75] (taken from Ref. [16]). 

advantage of the near-far detector approach to get a mo del- independent result for (V)SBL 
u e disappearance, because for a typical reactor neutrino energy of 1 MeV the oscillation 
length corresponding to Am 2 « 10 2 eV 2 is of the order of 1 cm. 

It is interesting to note that the near-far detector setup that we have chosen is sensitive to 
v e disappearance with small mixing (sin 2 29 > 2 x 10~ 3 ) for values of Am 2 as large as 10 3 eV. 
The condition for the observation of a spectral distortion caused by neutrino oscillations is 
that the uncertainty of the phase of the oscillations due to the energy resolution in Eq. (flOl) 
is smaller than about tt/2. One can easily find that this happens for neutrino energies 

-| 2/3 

(16) 



E > 



e Am 2 L, 



off 



2tt El 12 



where we have considered the effective baseline in Eq. ( Tl4|) . Since for the near detector 
L cS ~ 180 m, if Am 2 = 10 3 eV the condition ( Tl6l) is satisfied for E > 18 GeV. Since 
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Figure 6: Fits in the sin 2 29-Am 2 plane for three chosen test-points marked by the diamonds (la, 2a, 
3a, 2 d.o.f.). Here CPT invariance is assumed. Near (50 m) and far (2 km) detectors are used with our 
systematics model, the effects of averaging over the decay straights are taken into account. The unshaded 
contours show the result without averaging effects over the straights. They are too small to be visible in 
the middle panel. 

for the assumed = 25 GeV the neutrino energy spectrum extends up to 25 GeV, as 
shown by the curve in Fig.l of Ref. [67] with off-axis angle 9 = 0°, the oscillations are not 
completely averaged out in the highest-energy bins. This is illustrated in the left panel of 
Fig. HI in which the line corresponding to Am 2 = 10 3 eV has the constant averaged value 
0.5 sin 2 29 — 1 = —0.05 (for the assumed sin 2 29 = 0.1) only for E < 10 GeV. Other curves 
illustrate the distortion of the event rate spectrum for smaller values of Am 2 . One can see 
that the lower limit of the sensitivity to Am 2 of the near detector is about 1 eV , which 
instead produces a strong spectral distortion in the far detector (right panel of Fig. HJ). 

We also show in Fig. [5] a comparison with the default setup in Ref. [29] (thin dashed curve). 
This setup uses a low-gamma (7 ~ 30) beta beam using inverse beta decay as detection 
interaction, which means that it is not surprising that our result is about an order of 
magnitude better. Compared to Ref. [29], which uses only one detector and therefore runs 
in the systematics limitation in the larger Am 2 range, we also have very good sensitivity for 
large Am 2 . While both approaches rely on near detectors receiving neutrinos from a storage 
ring, they are conceptually very different: Ref. [29] uses the fact that the inverse beta decay 
reaction is well known to control systematics, whereas we control the shape error with two 
sets of detectors in the fashion of the new generation of reactor experiments. 

It is interesting to examine not only the sensitivity of our experimental setup to (V)SBL v e 
disappearance, which corresponds to a negative result producing an exclusion curve as that 
in Fig. [51 but also what could be the results if a signal is observed, i.e., v e and v e disappear. 

In Fig. [HI we show three qualitatively different possible results for the test values of the 
neutrino oscillation parameters marked by the diamonds. In the left panel, no degenerate 
solutions are present, and the parameters can be very well measured. There is hardly an 
effect of the averaging over the decay straights, as one can read off from the differences be- 
tween the shaded and unshaded contours, because the far detectors dominate the sensitivity 
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and oscillations have not yet developed at the near detectors. In the middle panel, we still 
have an excellent measurement dominated by the near detectors. In this case, however, the 
averaging effects over the straights are very important, and the contours without averaging 
are hardly visible. In particular, a degenerate solution appears at a smaller Am 2 . In the 
right panel, we show an even more extreme case, where only at the 2a confidence level 
sin 2 29 = can be excluded. 

5 CPT violation 

In this section we discuss the potentiality of the experimental setup described in Fig. [H 
with two pairs of near-far detectors, to reveal a violation of CPT symmetry, considering 
the different electron neutrino and antineutrino survival probabilities in Eqs. (j2J) and ([3]) as 
functions of the CPT asymmetries in Eq. (j5J). 

Since there are four independent parameters, given by Eqs. (EI) (J9]) , for simplicity we consider 
three test points inspired by Refs. [16-18]: 



Tl : 


sin 2 29 


= 0.05, 


Am 2 = 1.8 eV 2 , 


(17) 


T2 : 


sin 2 29 


= 0.1, 


Am 2 = 20 eV 2 , 


(18) 


T3 : 


sin 2 29 


= 0.1, 


Am 2 = 330 eV 2 , 


(19) 



and ocpt = ^cpt = 0. We fit the corresponding simulated data allowing for non-zero 
values of ocpt and tticpt in order to explore the sensitivity to the measurement of these 
parameters. 

The test point Tl is motivated by the best-fit of the data of the Bugey SBL reactor experi- 
ment [74] , which is compatible with the data of the Chooz reactor experiment [75] and the 
neutrino oscillation explanation of the Gallium anomaly [16]. The test points T2 and T3 are 
motivated by a possible explanation the MiniBooNE low-energy anomaly through VSBL v e 
disappearance which is compatible with the neutrino oscillation explanation of the Gallium 
anomaly [17,18]. Even if values of Am 2 larger than about leV 2 are incompatible with the 
existing standard cosmological bound on the sum of neutrino masses [9,76], we think that 
it is wise to test such bound in laboratory experiments. A violation of the bound may lead 
to a discovery of fundamental new physics related to non-standard cosmological effects. 

The best-fit regions for the three test points are shown in Fig. [71 The dashed curves represent 
the results without taking into account the averaging over the decay straights. The test-point 
Tl (upper row), with a relatively small Am 2 , is dominated by the far detectors, whereas 
in the near detectors (almost) no oscillations are present. Therefore, the cross sections can 
be directly reconstructed from the near detectors, and the fits are very clean. The effects 
of averaging over the straights are small because the signal sits in the far detectors, which 
sees a point source. The oscillation parameters can be measured at the level of 2% (la), 
and the CPT invariance can be constrained at the same level. 

The test-point T3 (lower row of Fig. [7J , is dominated by the short baseline, which means that 
the averaging effects over the straights are very important. The longer baseline measures 
the product of cross sections and 1 — 0.5 sin 2 29, which means that Am 2 can, before the 
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Figure 7: Best-fit regions in the Am 2 -sin 2 29 and acpT-^CPT planes for the test points defined in the 
main text (la, 2a, 3a, 2 d.o.f.). Note that different baselines where chosen for the test points. The dashed 
curves represent the results without taking into account the averaging over the decay straights. 
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Fi gure 8: Discovery reach for CPT violation from ocpt (left panel) or tocpt (right panel) as a function of 
the true sin 2 29 and true Am 2 . The different contours indicate for how small (true) values of acpT > (left) 
or mcpT > (right) CPT violation will discovered at the 3a confidence level, as labeled at the contours. 
The dashed curves show the result if the averaging over the decay straights is not taken into account. 

averaging over the straights (dashed curves), be very well measured compared to the mixing 
angle since it remains as a net effect between the two detectors. Only after the averaging 
effects, both oscillation parameters can be measured at the level of 1% (lcr), and the CPT 
invariance can be constrained at a similar level. 

The test-point T2 (middle row of Fig. [7]) shows a complicated case with an intricate inter- 
play between systematics and oscillation parameter correlations. Since there is an oscillation 
effect in both baselines, this case does not correspond to a classical near-far detector com- 
bination. The a priori excellent precisions for the oscillation parameters are spoilt by some 
complicated correlations. Nevertheless, percent level precisions are possible. 

Instead of constraining CPT invariance, we can also discuss the discovery reach for CPT 
violation. In this case, we assume that nature has implemented a small (positive) acpT 
or mcpT, and we fit the simulated data with the fixed parameters ocpt = "^cpt = 
(corresponding to CPT invariance), while we marginalize over the oscillation parameters 
sin 2 29 and Am 2 . We show in Fig. [SJthe discovery reach for CPT violation from ocpt (left 
panel) or mcpT (right panel) as a function of the true sin 2 29 and true Am 2 . The different 
contours indicate for how small (true) values of ocpt > (left) or m CPT > (right) CPT 
violation will discovered at the 3cr confidence level, as labeled at the contours. The dashed 
curves show the result if the averaging over the decay straights is not taken into account. 

From Fig. [HJ CPT violation may be discovered even if it is as small as 10~ 3 , provided that 
sin 2 29 is large enough. However, even for very small sin 2 29, a CPT violation of order 
unity is testable with our setup. Note that for larger Am 2 and especially for mcpT, the 
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Figure 9: Lower limit for |a C pT| obtained from Eq. (jSTJ with A^ T < -0.08, which is the 95% C.L. 
limit found in Ref. [18]. 



averaging over the decay straights strongly reduces the performance (by about one order of 
magnitude). 

In Ref. [18], a difference of A^ T = P ee - P ge - = — 0.17±g;g? at 90% C.L. was identified 
as the asymmetry between the electron neutrino and anti-neutrino VSBL disappearance 
probabilities which can explain the Gallium radioactive source experiments anomaly [14] 
and the MiniBooNE low-energy anomaly [13] without conflicting with the absence of P e 
disappearance in reactor neutrino experiments (see Ref. [77]). It is interesting to investigate 
if such CPT violation can be measured in a neutrino factory experiment with the near-far 
pairs of detectors that we have considered so far. 

Since in Ref. [18] Am 2 was considered to be large, in the range 



20 eV 2 < Am 2 < 330 eV 2 



(20) 



the neutrino and antineutrino survival probabilities were assumed to be averaged, leading 
to Ag e PT = 0.5 (sin 2 29p — sin 2 29 u ). In this case, the asymmetry ocpt is given by 



Ocpt = arcsm 



sin A9 j 



(21) 



Since IczcptI < 1, the mixing angle has a lower limit which depends on the value of A^ T . 
Moreover, since | sin4#| < 1 and 9 < tt/2, also IocptI has a lower limit, which is plotted in 
Fig. M for A^ T < -0.08, which is the 95% C.L. limit found in Ref. [18]. One can see that 
the bound on A^ T implies that sin 2 29 >4x 10~ 2 and IocptI ^ 0.10. Confronting these 
values with the left panel in Fig. [HI and taking into account that we consider the large values 
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of Am 2 in the range (|2DJ, it is clear that the CPT violation required by A^ T < -0.08 will 
be easily discovered in a neutrino factory experiment with the near-far pairs of detectors 
that we have considered. 

6 Summary and conclusions 

In this work we have discussed the potentiality of testing Short-BaseLine (SBL, with 10 _1 < 
Am| BL < 10 eV 2 ) and Very-Short-BaseLine (VSBL, with 10 < Am^ SBL < 10 3 eV 2 ) electron 
neutrino disappearance in a neutrino factory experiment, based on the current setup of 
the International Design Study for the Neutrino Factory (IDS-NF) [63]. Since this setup 
uses both muon and anti-muon decays, a possible difference between the neutrino and 
antineutrino disappearance can be studied, which could constitute a revolutionary discovery 
of CPT violation. 

We showed that for these purposes the ideal configuration would be the two pairs of near- 
far detectors (shown in Fig. [1]) in a similar fashion to reactor experiments with near and 
far detectors (Double Chooz [70], Daya Bay [71], etc.) to cancel systematics. The near 
detectors are chosen to be at a distance of about 50 m from the muon storage ring, in 
order to be sensitive to oscillations due to a Am 2 as large as about 10 3 eV 2 . For the far 
detectors an appropriate distance from the muon storage ring is about 2 km, which gives 
a good sensitivity to oscillations generated by a Am 2 as small as about 10 _1 eV 2 . In this 
way, it is possible to explore (V)SBL v e and v e disappearance with effective oscillation 
amplitude sin 2 29 as small as about 10 -3 for Am 2 > 1 eV 2 (see Fig. [5]) taking advantage of 
the comparison of the event rates measured in the near and far detectors, which reduces 
dramatically the systematic uncertainties due to insufficient knowledge of the cross sections 
(see the discussion in Section [3]). 

We have also shown, in Sectional that the chosen detector setup provides a good sensitivity 
to the measurement of a difference of the rates of v e and v e disappearance which would be 
a signal of CPT violation. For instance, our setup is sensitive to an asymmetry between the 
neutrino and antineutrino mass squared differences at the level of up to 10~ 3 , depending on 
the value of the mixing angle. Let us emphasize that a discovery of CPT violation would 
represent a revolution in our knowledge of fundamental physics, because the CPT symmetry 
is a fundamental symmetry of local relativistic Quantum Field Theory. Therefore, pursuing 
this line of investigation is of fundamental importance. 
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